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We analyze the time evolution of quantum entanglement in a model consisting of two two-level 
atoms interacting with a two-mode electromagnetic field for a variety of initial states. We study 
two different coupling schemes motivated by the forms that can arise due to atomic separation. We 
observe a variety of qualitative features such as entanglement sudden death, dynamical generation, 
protection, and transfer between subsystems. Our quantitative analysis shows that these cases with 
different couplings and initial states differ significantly in these qualitative features. The multifarious 
behaviors in these two-mode cases suggest the importance of considering atomic separation carefully 
for any model where two atoms interact with a common field. 



I. INTRODUCTION 

Quantum entanglement has been extensively studied, both due to its fundamental significance in quantum theory 
[U |2] and its utility as a resource for quantum communication and quantum information processing [31 H] . Atomic 
physics offers a domain with sufficient control of the system and isolation from noise that it has been possible to 
do extensive experimental study of quantum entanglement [SJ [B], and, therefore, it is also a productive target for 
theoretical study of the issue as well. 

One of the simplest scenarios for theoretically studying the dynamics of entanglement between atoms is that of two 
atoms which are isolated from one another and interact with different electromagnetic fields. Studying this sort of 
model led to the discovery of the entanglement sudden death (SD) phenomenon [711, in which entanglement decays 
to zero in a finite time rather than asymptotically. SD has garnered significant interest because it is both unintuitive 
and potentially undesirable. 

An alternative simple model for atom-field interaction in which to study entanglement dynamics is a Dicke model 
|12| .where one assumes all atoms are grouped in a sample whose size is small compared to the resonant wavelength 
(resulting in identical coupling to every atom) . However, such a simplistic model may miss the variety of behavior that 
can result when the atoms are not confined to such a small sample. Entanglement dynamics have been shown to have a 
significant distance dependence when two atoms are interacting with a common field [13 . For atoms weakly interacting 
with a continuum of field modes in the Born-Markov approximation, it has been shown |14fll7| that changing the 
atomic separation of two atoms can affect whether there is SD and whether there is revival of entanglement, as well 
as modify the dynamical generation of entanglement; in short, the qualitative features are sensitive to the atomic 
spacing. At short interatomic distance non-Markovian effects associated with induced interactions between the atoms 
due to the quantum field become more pronounced whose qualitative behavior varies greatly with different classes of 
initial states, as detailed in the study of [TH]. So it is clear that distance dependence gives rise to a significant variety 
of behaviors, and some seemingly innocuous approximations can qualitatively alter the entanglement dynamics in 
unintuitive ways [191 EO] . 

In view of this, we seek to study entanglement dynamics in a model that is sufficiently complex to manifest some of 
this variety of behavior yet simple enough that the dynamics may be understood in considerable detail and obtained 
with fewer approximations. The simplest model one may pursue along this line would be two two-level atoms (2LA) 
coupled to a single mode with couplings that reflect the distance dependence; however, as we demonstrate in Sec. [IT] 
one needs to include at least two field modes before any non-trivial distance dependence can arise in the problem. So 
we will study such a model with two field modes here. Out of the class of Hamiltonians that can arise from distance 
dependence, we will focus on two special cases with the aim of illustrating the variety of different behaviors that can 
arise. Specifically, we will show how SD, dynamical entanglement generation, and other phenomena differ between 
these two cases. We will also compare with the behavior of the analogous versions of the two well-studied types of 
models mentioned above: For two isolated atom-field systems, we will compare with the case where each system is a 
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single 2LA interacting with a single field mode |10l [TTl I21| . In the instance of two atoms interacting with the same 
field, we will compare with the properties of two 2LAs interacting with a single common field mode ^2^ [55] . We will 
show that the entanglement dynamics in these models differs significantly from the model of two atoms interacting 
with a two-mode field that is our focus. 

The remainder of this article is structured as follows: In Sec.|IT]we describe the models we use and the assumptions 
behind them. In Sec. |III| we describe the quantum states of interest. In Sec. |IV| we describe the cases where entangle- 
ment is dynamically generated from an initially separable state and where entanglement is transferred between the 
atom and the field sectors. In Sec.|V]we describe the cases where entanglement meets sudden death (SD) and identify 
the conditions under which entanglement can survive. Situations in which entanglement never dies (i.e., the system is 
never separable) we term always alive (AL). In addition to known scenarios of entanglement death, birth and revival 
we also encounter situations where entanglement dies only for an instant (DI). The qualitative features for all the 
cases we have studied are summarized in Table [Tj Finally in Sec. |VI| we conclude with a summary and discussion. 



II. THE MODEL 
A. System Hamiltonian 

Consider a pair of two-level atoms (2LA) coupled via the multipolar interaction Hamiltonian to a collection of 
electromagnetic (EM) field modes in the dipole and rotating- wave approximation [21]. The interaction Hamiltonian 
of the system in the Schrodinger picture has the form 

= ^J2Y1 dj.q^t^i + dlq^Jal, (1) 
j 1 

where and aj are the raising and lowering operators for the two-level system representing the j'^ atom, and 
and CLq are the creation and annihilation operators corresponding to the g*^ normal mode with frequency ujq in some 
mode decomposition of the EM field with classical electric field mode function (r). If we assume the atoms to 
have fixed center of mass positions Rj , assume that the atomic transition for each atom does not alter the angular 
momentum of the atom, and assume the dipoles associated with all the atomic transitions have identical direction e^, 
then 

= djBd ■ fg (Rj) y^^, (2) 

where the field is quantized in a region of volume V and dj is the complex dipole matrix element for the transition in 
the atom [25j. 

In general, the mode functions (Rj) can be quite complicated, and even for a single atom there can be position 
dependence in the dynamics arising from the boundary conditions that the mode functions obey. In order to distinguish 
these effects from the effect of atomic separation, we will consider traveling-wave mode functions that are solutions 
to an EM boundary value problem which is translation invariant in the coordinate separating the two atoms. For 
simplicity, we will discuss traveling plane waves in free space, though another example would be the modes of a 
toroidal resonator. 

Having selected traveling plane waves as our mode decomposition, we can then choose to define them such that for 
each wave vector k one of the polarizations is perpendicular to Bd, so that it does not couple to the problem and can 
be neglected. With this assumption we know that for all remaining modes 

In general each dj may have a different complex phase, so 

arg(5i,9) = arg(dj)-)-fe, -iij. (4) 
However, it turns out that, without loss of generality, one may study a much smaller set of Hamiltonians. 
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B. Equivalent Hamiltonians 

Given a particular set of complex phases for the coupling constants gj^q, one may make a trivial basis transformation 

C/b = ]Je-*«'*'/2-Qg<.ata,^ (5) 

I s 

which simply amounts to redefining the reference for the phases of the atoms by — > a^e^^^^ and for the fields by 
aq — > CLqe^'^'^'' . After this transformation, the interaction Hamiltonian becomes 

H'l - U.HiUl =Hi=hY,Y. + idj;) * ^7^1 (6) 

with 

3;,,=5,,,e-«^e-<', (7) 
and we may obtain the solution to the original problem by using the transformed Hamiltonian: 

1^ (t)) = e~'"''/^ \ij (0)) = Ule-'"''''''Ub \^ (0)) = ule-'"'''''' (0)) . (8) 

We may always choose a set of phase shifts ^; and Cs such that for some particular j and m 

arg = arg(d,) + fes • -R; - 6 - Cs = if / = j or s = m (9) 

making those couplings real, but for other I and s the phase difference 

arg {gi^s) - arg {gj^q) = arg (d/) - arg {dj) + ks ■ Ri - kq ■ Rj (10) 

cannot generally be eliminated by the added phase shifts and represents the origin of non-trivial physical distance 
dependence. 

If we consider two identical atoms, with frequencies ujq and dipole strengths \dj\ — d coupled to a single EM field 
mode, then the previous paragraph implies that it suffices to consider only a Hamiltonian where both couplings are 
real, and the atomic separation does not enter; thus, there can be no non-trivial distance dependence. However, if we 
consider the two identical atoms coupled to two field modes, then without loss of generality we can write the total 
Hamiltonian of the system as 

H = Hljq ^(t|(Ti + ^2^2^ + hujid\di + huj2a\a2 + hgi (^a^di + (Tj~a| 



^.92 {(^t^2 + O"! 02 ) + %1 (o'^ai + (72 d\) + %2e*'^ ['^2^2 + 0'2 ^2 ) ' (H) 



where all distance dependence arises from (j) = (^2 — fci) ■ {R2 — Ri)- For simplicity, we will further assume that 
(jji ^ UJ2 = 1^0 and |fci • e^l ~ \k2 ■ e^j, which by Eq. (|2| implies gi — g2 = g. Our aim is to get a sense of the variety 
of different entanglement dynamics that can result from different separations and initial conditions. To that end, we 
will consider two special cases, which are arguably extreme cases of the general model, ^ = and (f> — tt. 

C. Mapping Equivalent Models and Time Evolution 

Since we have assumed the two field modes have the same frequency, rather than using the original modes Fi and 
F2 one could equally well choose a different mode decomposition of the field where the two modes are replaced by 
linear combinations TFi and TF2 with annihilation operators 

11 = -^(ai+a2) (12) 

12 = ^ (ai - 02) . (13) 
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In the case of two modes with symmetrical coupling (TMSC), where cj) = 0, the interaction Hamiltonian can then be 
written 



(14) 



in terms of these transformed modes, so that the atoms only couple to Ai and not ^2- This shows that the Hamiltonian 
is equivalent to the model of a single mode symmetrically coupled (SMSC) to two atoms. If we consider the evolution 
of the reduced density matrix of the atom pA (t) in the TMSC model it should be the same as in the SMSC model with 
a properly transformed initial state. Namely, if the total system (atoms and modes) is in an initial state described by 
the density matrix xtmsc (0) j then the appropriate initial density matrix for the SMSC model is obtained by making 
the mode transformation of Eqs. (12 1 and (13 1 and tracing out the second transformed mode Ti^2- For example, if 



the initial state in the TMSC model is separable with the field in a product of Glauber coherent states |a,/3), then 

,,. a + /3 a-/3\/a + /3 a - P 



XTMSC (0) = 10) 



|a,/3) {a,f3\ Tr 



TFn 



/2 ' V2 

a + f3\ /a + f3 



V2 



V2 



XSMSC (0) (15) 



is the appropriate mapping to the equivalent SMSC problem. It is important to note that, because this mapping of 
initial states involves a partial trace, it is a many-to-one mapping from the TMSC problem to the SMSC problem 
(and this mapping does not preserve purity). In order to solve the dynamics in the SMSC model, and by extension 
the TMSC model, we simply compute the time evolution operator expressed in the atomic basis {|ee) , \eg) , \ge) , \gg)} 
directly by exponentiation (as in, e.g. |23|V 



Hi = V2g 
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where 



Si = Ai 



sin (^\/AAgt^ 
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(16) 



When (f) — TT we have a two- mode model with asymmetrical coupling (TMAC), and we may again use the mode 
transformation of Eqs. (12 I and (13 1 and write the Hamiltonian as 



Hi = V2hg a^Ai + &^ A\ + A2 + A\ 



(17) 



In this case, rather than both atoms coupling to one mode we see that atom one couples only to transformed mode 
TFi while atom two couples only to TF2 , implying that this Hamiltonian is equivalent to a model comprised of two 
subsystems that are totally isolated from one another, each composed of a single atom coupled to a single mode. We 
will call this the double Jaynes-Cummings (DJC) model. This sort of model with isolated subsystems is common to 
the study of entanglement sudden death [3 [5], and the DJC model specifically has been studied pUl ITTl [2T| . As in 
the previous case, the evolution of pA (t) for the TMAC model should be the same as given by the DJC model with 
the proper mapping of initial states. In this case the mapping of initial states is limited to transforming the modes 



according to Eqs. (12) and (13 I, so 



|a,/3) 
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This mapping implies that there can, for example, be no dynamical generation of atomic entanglement in the TMAC 
model unless the DJC initial field state obtained by the mapping is entangled. In the DJC model we can write the 
unitary time-evolution operators for the two separate non-interacting atom-field subsystems as Ui and U2, and then 
the total time evolution operator is U = Ui ® U2- We again compute the two subsystem unitary time evolution 
operators by direct exponentiation to obtain 
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D. Quantifying Entanglement 

In order to quantitatively study the entanglement dynamics, we must specify a way of quantifying entanglement. 
For two qubits one may compute the entanglement measure known as the entanglement of formation efficiently in 
terms of the Wootters concurrence f^. In order to compute the concurrence, one must choose an arbitrary basis and 
construct the spin-fiipped operator 

p={ay®(jy)p*{ay®ay), {22) 

where the star denotes the complex conjugation in that basis and Gy is the Pauli matrix in the same basis. The 
concurrence is 

C(/5) =max|0, v/A^- v/A^- v/a;^- V/A4}, (23) 

where the AjS are the eigenvalues of the matrix pp in decreasing order. The entanglement of formation can then be 
computed as 

^HC)-k(^-±^), (24) 

where 

h {x) = -X log2 [x) - (1 - x) log2 (1 - a:) . (25) 

Alternatively, to quantify entanglement one can also use a monotone known as the negativity [27' (which quantifies 
the degree to which the state violates the Peres-Horodecki positive partial transpose condition), defined 

^{p)=^^^ , (26) 

where jf"^ is the partial transpose of the density matrix, defined by 

{j,k\F'\i.q)^{j.q\p\i.k) (27) 

in an arbitrary basis B, and ||-||-^ is the trace norm. Either of these approaches can be used to quantify entanglement. 
The negativity is arguably less precise, because in some systems it can be zero even for an entangled state, but we 
will not consider such situation. The advantage of the negativity is that it can be efficiently computed even in a large 
Hilbert space. We will most often use the concurrence in this article; however, we do quantify the entanglement of 
field states using the negativity. 
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III. THE STATES 



In order to illustrate the variety of behavior that can arise among the four models we discussed in Sec. II C we 
will examine the entanglement dynamics of a selection of initial states in which the atoms are separable from the 
fields and which are comprised of familiar atomic and field states. In each case we will select the atomic state from 
the set of pure states {\gg) , |ee) , \eg) , |4>) , l^-)}, where |$) = (|ee) + \gg)) /V^, l^-) = {\eg) + \ge)) /v^, and e and g 
label the excited and ground states of the atom respectively. The initial field state will be either the vacuum |00), a 
product of Fock states |njv, m^v), a product of Glauber coherent states \ac, /3c)i a product of squeezed vacuum states 
\S,sq, ~^sq) , a two-mode squeezed vacuum state (TMSS) 0, 0) , a thermal state pth (with both modes having equal 
temperature), the pure state \r]nm)i or the mixed state Pnm- The two-mode squeezed state is defined as the state 
resulting from the action of the two-mode squeezing operator S'(^) = 0102-40^02) vacuum. The state |?7„m) 
is the result of mapping the state |njv,mAr) in the original modes of the TMAC problem to the transformed modes 
equivalent to the DJC problem, with 

n m 

and Pnm = TrTi?2 [IVnm) iVnmW IS the state in the SMSC model that gives equivalent evolution to the state In^r, mjv) 
in the TMSC model. "Cfe represents the binomial coefficients. 

The correspondence drawn between the TMSC-SMSC leads to the essential feature that the map for the initial 
field states from TMSC to SMSC is many to one; the set of initial states that have the same reduced density matrix 
for the first transformed mode TFi have identical entanglement dynamics in terms of A1-A2 entanglement. As a 
counterintuitive example of this feature we will see that a squeezed state of the form \^sq, ^S.sq, 0) and a thermal field 
in the TMSC model give the same entanglement dynamics provided the average number of photons in the thermal 
field corresponds to that in the squeezed state {fith ~ sinh^ l^sijD- 

We summarize our findings for the entanglement behavior given the various initial states considered in the four 
models in Table [T) listing the equivalent TMSC-SMSC cases and TMAC-DJC cases. When discussing entanglement 
sudden death, we adopt the usage of the term as in [IT in applying it only to instances where the entanglement goes 
to zero for some time interval of non-zero length. In the case where entanglement goes to zero only for an instant 
during the time evolution we refer to it as death for an instant (DI). If there is a non-zero entanglement at all times 
once it is generated in the system then we label it as being "always living" (AL). 



IV. ENTANGLEMENT GENERATION AND TRANSFER 



The generation of entanglement from an initially separable state by the dynamics of the system is one interesting, 
and potentially useful, phenomenon to examine. In the DJC model each atom interacts only with a separate field 
mode, so the dynamics cannot increase entanglement between the two atom-field subsystems; therefore, if the atomic 
state is not entangled initially then it will remain separable, unless there is an initial entanglement between the field 
modes that can be transfered to the atoms by the dynamics. Knowing this, we can see that any initial field state for 
the TMAC model that maps to a separable DJC field state will also fail to generate entanglement. The nature of the 
mapping means that even some entangled field states will fail to generate atomic entanglement in the TMAC model, 
while some separable states will map to an entangled DJC state and will generate entanglement. 

This structure means that, as shown in Table |T] many familiar initial field states fail to dynamically generate 
entanglement in the TMAC model. In contrast, entanglement generation is a common feature in the TMSC model 
for the same selection of field states. For example, starting with a two-mode squeezed state in the field modes does 
not generate any entanglement in TMAC, because the corresponding field state in DJC possesses no correlations 
amongst TFi and TF2. On the other hand, a two- mode squeezed state in the symmetrically coupled case generates 
entanglement in the initially separable atomic states as can be seen in Fig. [l] In fact, for initial atomic states |ee) 
and \gg), if the field is sufficiently squeezed then entanglement once generated sustains forever. 

When considering the TMAC model with an initial field state that is a product of squeezed states \£^sq, —^sq) we 
find that entanglement is dynamically generated as shown in Fig. [2] One can understand this from the fact that 
the state j^^g, —£,sq) maps to the DJC model with an initial TMSS, so the generation of atomic entanglement occurs 
simply because the dynamics transfers the entanglement between the field modes to the atoms. Further, we find that 
there is an optimal squeezing value that generates maximal entanglement. The maxima in generated entanglement 
occur for values of the squeezing parameter such that the field state is close to being a maximally entangled qubit 
state (ao(^) |00) +ai{^) |11)), with ai(^) and ao(C) being comparable. On increasing the squeezing parameter further 
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SD 
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No 


No 


SD 


SD 


6. \isq, -isq) 
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SD 
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No 
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^ No entanglement for rijv = "m-N 



Table I: Entanglement dynamics for two modes symmetrical coupling (TMSC) with = and anti-symmetrical 
coupling with (f> = tt. Columns A-C list whether there is an entanglement generation in an initially separable atomic 
state (yes or no) . The dynamical phenomena observed in columns A-E are listed as entanglement sudden death 
(SD), entanglement dies for only an instant (DI), entanglement remains non-zero at all times and so is "always 
living" (AL). A '/' denotes that both kinds of dynamics are present depending on the particular initial state chosen 
from within the class indicated. The initial states have been explained in Sec. |III 



there are contributions from higher Fock states which decrease the transfer of entanglement. We also notice some 
secondary peaks on increasing the squeezing parameter ^^q- 

By comparing the effects of a TMSS and thermal field state in the DJC model, one can also get some insight into 
the role of correlations between the field modes. The reduced density matrices for individual modes correspond to a 
thermal state for both the cases; the difference being that in a TMSS the two field modes are strongly correlated with 
each other while in the thermal fields there are no correlations. One would generally expect that the correlated TMSS 
generates more entanglement in the atomic subsystem, which is trivially true for the case of an initially separable 
atomic state. This intuition generally extends to the situation of initially entangled atomic states, where we see that 
apart from the regular SD pattern in the absence of the field-field correlations there are spontaneously generated 
peaks as a result of two-mode squeezing(Fig. [S]); however, we note as an exception that at certain instants of time the 
entanglement in the presence of an initial thermal field exceeds that of the TMSS. Transfer of field-field {TF1-TF2) 
correlations can be explicitly seen in a small squeezing approximation when a two mode squeezed state interacts with 
an initially separable atomic state |ee) (Fig.|4|. Since the probability of higher photon numbers is small, by restricting 
to the 4-dimensional subspace of lowest energy states in the Fock basis for TF1-TF2 we obtain the negativities for the 
two subsystems as, 

Na.^a, ~\ min(5?c? - ^s\clM (29) 
AfTFi~TF2 «|min(siCi - £,clcl,0)\. 

where si = sm{^/2gt), ci = cos{\/2gt) and C2 = cos{2gt) 

In the TMSC model, the same initial field state, consisting of a product of squeezed states, corresponds to having 
a thermal field state in the SMSC model, a situation where entanglement generation has been previously studied 
|23| . Comparing the entanglement behavior in Fig. [5] for the two different couplings, we see that for the initial 
atomic state \gg) there is more entanglement generated in TMAC than TMSC. As mentioned before, two thermal 
fields in TMSC also map to the same situation. This is a somewhat unintuitive situation for the generation of 



entanglement. There is no entanglement generation when the atoms are initially in the state |ee). As Figs. 6c and 6d 
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Figure 1: Entanglement dynamics for a single mode squeezed state (|^s)) in SMSC or two-mode squeezed state 
(1^,0,0)) in TMSC interacting with different initial atomic states 



show, for \eg) entanglement is generated with only instantaneous zeros, however high the temperature, and for \gg) 
there is entanglement generation and sudden death behavior with an optimal fith ~ 1 for the generation of atomic 
entanglement. 

The observation that the initial atomic state |ee) will remain separable can, in fact, be generalized to any initial 
field state with a density matrix that is diagonal in the Fock basis. As had been shown in [5^, in the SMSC model 
an atomic state |ee) interacting with any Fock state \n) in the field mode never gets entangled This then implies 
that the atom-atom density matrix given as p^"^^ = Trj? [IV'n) (V'nl] remains separable, where the time evolved state 
I ■(/'„) = U \ ee) \n). Extending to a general completely- mixed field density matrix the time evolved atom-atom density 
matrix given as pi2 = Tr^ \J2n-^n-P''"'^\ clearly a convex sum of separable density matrices, and hence there is no 
entanglement generation. 



Looking at the initial state \eg) \n) in the SMSC case, we observe no SD in Fig. 6a This can be explained by 
considering the state as a superposition \eg) — ^ {\eg) + \ge)) ~\- | {\eg) — |<?e)), where due to the symmetry of the 
coupling constants the maximally entangled dark state {\eg) — \ge)) /\/2 does not interact with the field. It is only 



To explain this it can be observed from symmetry arguments that the time evolved atom-atom density matrix is an incoherent mixture 
of the states {|ee) , \gg) , |^)} such that the contribution of the maximally entangled part |^) in the mixture is at all times smaller as 
compared with the other two. 
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(a) |ee) - Entanglement generation from transfer of (b) \eg) - Secondary peaks in entanglement vs squeezing can 

correlations from the field, optimal squeezing can be seen be observed for ^ 1.5 

^^1 



Figure 2: Entanglement dynamics for a two-mode squeezed state (|^,0,0)) in DJC or product of single mode 
squeezed states (|<^s, —£,s,0)) in TMAC interacting with different initial atomic states 



TMSSin DJC, ^=0.5 
Thermal states in DJC, njj^=sinh''(^) 




Figure 3: Comparing entanglement dynamics for a two-mode squeezed state vs a thermal state in the DJC model for 

an initial atomic state |$) 



momentarily during the evolution that the state of the system returns to being the original separable superposition 
\eg). As a result we always observe some entanglement between the two atoms for an initial field density matrix 
diagonal in the Fock basis. The dynamics of the \gg) \n) state in the SMSC model, shown in Fig. 6b exhibit SD in 
general except for the special case of n = 1 where because of symmetry reasons the state oscillates between the states 
\gg) |1) and \'^) |0), going from being separable to maximally entangled. Hence, any density matrix diagonal in the 
Fock basis with a high component of |1) (1| would generate more entanglement in general. 



Apart from the Fock state in the SMSC model with a thermal field, another example of having a initial field density 
matrix diagonal in the Fock basis is to have a Fock state \n, m) in the TMSC which corresponds to the SMSC density 
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Figure 4: Small squeezing entanglement dynamics for a two-mode squeezed state in the DJC model interacting with 
an initially separable atomic state |ee), entanglement is generated in the atomic subsystem via transfer of field-field 
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matrix 

^ n m 

Pn^n = 2m+n^|^| Z! \m + Jl ~ k - I) {m + n - k - l\ (30) 

fc,p=0 l,q=0 

where Kmnfc/ = "'Ck"Cp"^Ci"^CqSi^^ip+q{m + n — k — l)\{k + /)!(— 1)'. For this state again we see no entanglement 
generation for |ee) and DI for \eg) and maximal entanglement in \gg) |10), as shown in Figs, 
point we observe is that for an initial state \gg) |n, n) there is no entanglement generation, which is a common feature 
between TMSC and TMAC. A Fock state \n, m) in TMAC transforms into an entangled state in the DJC model, so we 
expect and observe entanglement generation in the system for an initially separable atomic state. As an exception we 
find that for n = to, if there is no atom-atom entanglement to begin with then the atoms remain separable evolving 
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into a completely mixed state. This is counterintuitive in the sense that even though the field state in the DJC 
model is entangled to begin with, there is still no transfer to the atomic subsystem. This feature can be explained by 
considering the initial field state |77„„) = ^^^^ 12k=o "Cfe(^l)'^\/2/c!(2n — 2k)\ \2n — 2k, 2k). If the atoms are initially 
in the state |ee) then time evolution will lead to an entanglement of atomic and field states such that detecting whether 
the number of photons in each mode is even or odd tells us the state of the two atoms. On tracing out the field this 
gives us a completely mixed atom-atom density matrix with no atom-atom entanglement. The same is true for the 
initial atomic states \eg) and \gg). 

While there are more nuanced details in what we have reported in this section, generically speaking it seems that 
the TMSC (and SMSC) models are much more conducive to the dynamical generation of entanglement than the 
TMAC (and DJC) model, suggesting that atomic separation may have a strong infiuence on this. Particularly useful 
is the ability in the TMSC to dynamically generate entanglement with a thermal field state or from a squeezed state 
in SMSC where one can even generate atomic entanglement that is AL. 



V. ENTANGLEMENT SUDDEN DEATH AND PROTECTION 



The phenomenon of entanglement sudden death has clearly provoked much theoretical interest, and it is related to 
another question that is both interesting from a theoretical perspective and clearly of great practical importance: how 
can one protect a system from disentanglement? Here we do not propose any active scheme for protecting entanglement 
(as in, e.g. [IH]), but rather consider what initial states of the field tend to minimize the loss of entanglement or 
safeguard entanglement once it has been dynamically generated. Of particular interest is avoiding SD. 

In terms of the effect of the initial atomic state on the entanglement dynamics, it has been discussed previously 
in [TD] for the DJC model that the initial state |<i>) |00) undergoes sudden death while I^P) |00) has DI of the atomic 
entanglement. This differentiation in behavior is common to many of the models for studying SD in which each atom 
interacts with a separate field [7, 8J. The same situation must occur in the TMAC case as well. In comparing these two 
initial states for the TMSC we find a sort of reversal of roles; as Fig. [t] shows, |$) |00) enjoys non-zero entanglement 
at all times and |\[') |00) still suffers DI, so now the entanglement of |$) is better preserved. This reversal also holds 
in the case of a thermal field in the TMSC model so long as the thermal average photon number is below a threshold 
value ficrit ~ 0.43 (Fig.[8| with the state being AL (above this critical temperature |$) experiences SD as well). 

If one considers the TMSC model with the field modes in a two-mode squeezed state, then for the separable initial 
atomic states |ee) and \gg); if the field is sufficiently squeezed, entanglement is dynamically generated and once 
generated sustains forever (whereas a Fock state or thermal state may generate entanglement but it goes to zero again 



at some later time). This behavior is shown in Figs, la and Ic If the atoms are initially in the entangled atomic 
state 1$), then Fig. |9] shows that in the TMSC model increasing squeezing raises the minimum value of entanglement 
progressively towards a situation where the state is maximally entangled at all times, while for the TMAC model the 
system exhibits SD, and the entanglement gets destroyed in general with increasing squeezing. 

A general observation that seems to fit for most of the cases considered is that a higher average number of photons 
in the field destroys entanglement. As an exception, however, it can be seen from Fig. |10a| that in the TMSC model 
with the field in a product of coherent states, if the average number of photons is in a particular range then there is 
no SD once entanglement is generated for the states |ee) and \gg). While Fig. |10b| shows that for \eg), the regular 
rule applies. In the case of the initial field being a two-mode squeezed state (Fig. Ill), we observe that after a threshold 
squeezing (or average photon number) there is no sudden death of entanglement after generation in the system when 
the atoms are initially in the state |ee) or \gg). On the other hand, for the initial atomic state \eg) entanglement 
decreases as the average number of photons is increased. 

Much as for generating entanglement, we find that the TMSC model is generically better suited to protecting 
entanglement from sudden death. While in the TMAC model at nonzero temperature disentanglement occurs for 
both entangled states considered, in the TMSC we find that below a threshold temperature entanglement of the 1$) 
state remains AL. In addition, we find that in the SMSC model a two-mode squeezed vacuum field can keep |$) 
almost maximally entangled at all times, provided squeezing is large enough. In the TMAC sudden death is a quite 
generic feature, which is only escaped for an initially entangled state when the field is the vacuum or a select product 
Fock state. 



VI. SUMMARY AND DISCUSSION 



We have analyzed the entanglement dynamics in a model consisting of two two-level atoms and two electromagnetic 
field modes for a variety of familiar field states and classified the various cases in terms of phenomena such as 
dynamical entanglement generation, entanglement sudden death. One aim of this analysis is to get a sense of the 
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variety of different classes of entanglement dynamics that can arise from different atomic separations in a case where 
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Figure 7: Entanglement dynamics in vacuum interacting with initially entangled atomic states |3>) and in 

presence of symmetric and anti-symmetric couplings 




Figure 8: Entanglement dynamics for a single mode thermal field interacting with an initially entangled atomic state 
1$) in the TMSC model. SD occurs after a certain threshold temperature. 



two atoms interact with a shared EM field. It is useful to examine this question in a very simple model that can 
be solved with a minimum of approximations (using which can lead to unphysical effects) and understood in detail. 
We have argued that there is no non-trivial distance dependence in a single mode model, and. therefore, a two-mode 
model represents the simplest case for the study of distance dependence in the entanglement dynamics. 

We have studied two arguably extreme cases out of the class of two-mode Hamiltonians that can arise, one in which 
the two modes are symmetrically coupled (TMSC) to the second atom, and one where the two modes have asymmetric 
coupling (TMAC). A useful insight in understanding these models is that the atomic dynamics in the TMSC model 
correspond exactly to the dynamics for a model with single field mode symmetrically coupled (SMSC) to both atoms 
with a suitable mapping of the field state, while the atomic dynamics for the TMAC model correspond exactly to 
the dynamics for a double Jaynes-Cummings (DJC) model, made up of two isolated subsystems each with one atom 
coupled to one mode, under the proper transformation of the field state. These mappings help one understand the 
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Figure 9: Comparing entanglement dynamics for the initially entangled atomic states interacting with a two-mode 

squeezed state and in TMSC and TMAC 




Figure 10: Entanglement generation in presence of a coherent state in SMSC 



significant differences in behavior in the two seemingly similar models, giving a window into how significantly atomic 
separation can affect entanglement dynamics. Another significant implication of the mapping between the TMSC 
and SMSC models comes from the fact that the mapping of initial field states between those models is many-to-one 
(because it involves a partial trace); this shows us entire classes of field states for the TMSC model that will give 
exactly the same atomic dynamics. In particular, we saw that a product squeezed state can be identical to a thermal 
state with respect to the atomic dynamics. 

In examining the dynamical generation of entanglement from an initially separable atomic state in Sec. |IV| we 
find quite a marked contrast between the TMSC and TMAC models. While entanglement generation is a relatively 
common feature, present for a variety of field states, in the TMSC model, it is comparatively much more rare in the 
TMAC model. One aspect that highlights these differences is that in the TMAC model entanglement can be generated 
by a product of squeezed states or fock states, but in the TMSC model it can also be generated by more easily prepared 
states including thermal states and product coherent states. This difference is not so surprising however, if one views 
it in terms of the mappings we have introduced to the other models. When considering that the TMAC maps to 
the DJC model, where the two subsystems are isolated, one would expect entanglement generation to be relatively 
rare; it can only exist in cases where the field state is mapped to an entangled state whose entanglement can then 
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be transferred to the atoms. With the TMSC, by contrast, we have a mapping to the SMSC model where a single 
shared field mode can readily introduce entanglement between the two atoms. 

When one is concerned with protecting the entanglement of two initially entangled atoms, our analysis in Sec. [V| 
shows that again the TMSC model is better for that purpose in most cases. In the TMAC entanglement sudden death 
(SD) is a fairly generic feature, with the initial state \'^) |00) being one of only two classes we consider that does not 
show SD; the |00) shows "death for an instant" (DI), where entanglement goes to zero only on a finite set of points. 
The TMSC shows a reversal of the fortunes of the |00) and |$) |00) states, with the former experiencing DI while 
entanglement for the latter is AL, staying non-zero for all times. Moreover, the DI property of the dynamics of the 
1 4") 1 00) state in the TMAC model is fragile, in the sense that it is destroyed by even the smallest departure from the 
vacuum to, for example, a finite temperature field, while for the TMSC model the AL feature of the dynamics of the 
1$) |00) state is robust, remaining for non-zero temperature below a threshold ficrit — 0.43. This gives us a condition 
for protecting entanglement in this case. 

A different sort of issue we have touched on briefiy is the role of quantum correlations between field modes in the 
entanglement dynamics of the atoms. In the DJC model we compared the resulting entanglement dynamics for fields 
in a two- mode squeezed state (TMSS) and a thermal state. For appropriate choice of the squeezing parameter, the 
TMSS has the same reduced density matrix for either mode alone as the thermal state, so in this sense these form 
a pair of extreme cases to compare, one with strong quantum correlations while the other has none at all. It is 
necessarily true in the DJC that a thermal state cannot generate entanglement, while a TMSS does. Intuition would 
also suggest that the local entropy destroys entanglement and leads to qualitative sudden death features, while the 
entanglement generation we observe in the case of a two-mode squeezed state arises from the field-field entanglement 
being transferred to the atomic subsystem. The transfer of entanglement from field-field to atoms was analytically 
verified for a small squeezing approximation of the two-mode squeezed state where one can observe the field-field 
entanglement going to zero as the atom-atom entanglement builds up. For initial entangled states we would expect 
that the correlations of the TMSS aid in maintaining the entanglement of an initially entangled state, and we find 
this to generally be the case, although in some cases for brief periods of time a thermal state can actually result in 
greater atomic entanglement than the corresponding squeezed states. 

As a general result for the class of density matrices that are completely mixed in the Fock basis in SMSC case, we 
conclude that in terms of generating entanglement it is preferable to choose an initial atomic state \eg), which has 
entanglement being AL or at the least DI, as opposed to |ee) where no entanglement is generated. This had been 
previously pointed out in the case of Fock states and thermal states [221 123 • 

As an interesting result in terms of entanglement protection, we find that a single mode squeezed state interacting 
with symmetrically coupled atoms initially in the state |<i>) can be extremely effective in protecting the entanglement. 
Even for the vacuum, this entangled state is AL, but as squeezing is increased the minimum entanglement rises 
monotonically toward maximal entanglement. For entanglement generation by a single mode squeezed state in the 
separable atomic initial states |ee) and \gg), we found that for a sufficient amount of squeezing entanglement is not 
only generated but remains AL for all future times. Similarly, for a single mode coherent state there is a range of 
values of the average photon number for which there the generated atom-atom entanglement is AL. In these cases, it 
was observed that \eg) shows sudden death as squeezing or the average photon numbers were increased. 

We have chosen in this analysis to try to isolate the effect of atomic separation on entanglement dynamics from 
the position-dependent effects arising from boundary conditions. For this purpose, we have assumed in Sec. |II A| 
an atom-field coupling that depends on the coordinate separating the atoms only by a phase factor. For clarity we 
have supposed our two field modes are traveling- wave modes in free space, however a more experimentally relevant 
situation would be a toroidal resonator, where there is rotational invariance in the azimuthal direction that satisfies 
our requirements. In this case the two modes of interest would be two resonant, counter-propagating whispering- 
gallery modes. Because strong coupling between an atom and a whispering-gallery modes of a microtoroidal resonator 
has observed experimentally [29], there is the possibility of experimentally probing quite directly the model we have 
considered. However, in the experimental system there will be dissipative dynamics arising from emission into other 
modes outside the resonator, so a detailed comparison would require either including the dissipative effects in the 
theoretical model or a restriction to the case of sufficiently strong coupling to the resonator modes and early times 
that the dissipation could be neglected. 

Our analysis of special cases of the entanglement dynamics arising in two atoms interacting with two modes suggests 
a wide variety of different behaviors can arise, with qualitative features of the dynamics changing entirely between 
the two cases, even with the same initial field state. This suggests that understanding the distance dependence of 
entanglement dynamics for multiple atoms interacting with a common field will be quite important for predicting 
even the qualitative features that may arise. Furthermore, if one has the practical goal of dynamically generating 
entanglement or protecting entanglement once generated, the special cases we have considered suggest that the ability 
to achieve these goals will be greatly impacted by the separation of the atoms. Having gotten some sense of the 
variation in behavior that can arise, the obvious next step would be to quantify the entanglement dynamics over a 
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significant range of atomic separations rather than just special or extreme cases. The variety of qualitative features 
of entanglement dynamics we have illustrated here argue for the importance of these considerations. 
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